Interaction induced delocalisation for two particles in a periodic potential 



O 

o 
o 

(N 

> 

O 

ON 



Julien Vidal 1 , Bcnoit Doucot 2 , Remy Mosseri 3 and Patrick Butaud 4 
1 Laboratoire de Physique des Solides, CNRS UMR 8502, Universite Paris Sud, Bdtiment 510, 91405 Orsay, France 
2 Laboratoire de Physique Theorique et Hautes Energies, CNRS UMR 7589, Universites Paris 6 et 7, 

4, place Jussieu, 75252 Paris Cedex 05 France 
3 Groupe de Physique des Solides, CNRS UMR 7588, Universites Paris 6 et 7, 
2, place Jussieu, 75251 Paris Cedex 05 France 
4 Centre de Recherches sur les Tres Basses Temperatures-CNRS associe a I'Universite Joseph Fourier 

25, avenue des Martyrs, 38042 Grenoble, France 

We consider two interacting particles evolving in a one-dimensional periodic structure embedded in 
a magnetic field. We show that the strong localization induced by the magnetic field for particular 
values of the flux per unit cell is destroyed as soon as the particles interact. We study the spectral 
and the dynamical aspects of this transition. 

PACS numbers: 7f.fO.Fd, 7f.23.An, 73.20.Jc 



i 

CO 



I 

o 
o 



(N 
> 

in 
(N 

in 
o 
o 
o 



c 

o 
o 



X 



As shown by Anderson in 1958 tl , a quantum parti- 
cle in a disordered potential may be trapped in spatially 
localized eigenstates. A natural question which arises 
is whether such a dramatic effect operates in a many- 
body interacting system. This issue has been strong! 
revived by a series of experiments on Si MOSFETs 
which seem to indicate a metal-insulator transition in 
two dimensions. Since then, a lot of experimental and 
theoretical activity _has been dedicated to this highly 
controversial topic . Given the complexity of the full 
many-body systems, various groups have-siudied the al- 
ready non trivial two-particle problem ta ae . In partic- 
ular, D. L. Shepelyansky has shown convincingly that 
some two-particle eigenstates exhibit much larger local- 
ization length than single particle ones. Note however 
that in other situations such as in quasiperiodic systems, 
interactions .may generate strongly localized two-particle 
eigenstates 0. In addition, it seems that for repulsive 
interactions these interesting effects do not appear in the 
vicinity of the two-particle ground state, and this leaves 
open questions for the many electron case. 

Recently, it has been shown that an extreme local- 
ization mechanism induced by the magnetic field can 
lead to a complete confinement of the particle motion 
inside Aharonov-Bohm cages El . Contrary to the Ander- 
son localization, this phenomenon occurs in pure two- 
dimensional systems, i.e. without disorder, and is due 
to a subtle interplay between the structure geometry 
and the magnetic field. Two series of experiments have 
confirmed the existence of these Aharonov-Bohm cages. 
In the first one o , superconducting wire networks with 
the adapted structure exhibit a striking reduction of 
the critical current for the predicted varies of the mag- 
netic field. The second experiments E3 measure the 
magneto-resistance oscillations in two-dimensional meso- 
scopic structures with a small number of conduction 
channels and large electronic mean free path. A clearcut 
dip at half a flux quantum per loop is observed and con- 



firms the presence of this peculiar localization process. 

It is therefore very natural to ask whether these cages 
survive for interacting particles. In this letter, we present 
an exacly solvable model for two interacting particles un- 
der a magnetic field. To simplify, we deal with a quasi- 
one-dimensional model which exhibits Aharonov-Bohm 
cages. We show that for half a flux quantum per loop, 
dispersive two-particle bound states appear even for re- 
pulsive local interaction. In this system, the two-particle 
ground state is non dispersive but the first dispersive 
band is rather close in energy. Slightly away from these 
remarkable fluxes, these bound states survive until they 
merge in a two-particle continuum. Finally, we are led to 
speculate that a finite repulsive local interaction is able 
to turn the fully localized non interacting system into a 
strongly correlated metal provided the electron density 
is large enough. 

We consider a one-dimensional chain of square loops 
with periodic boundary conditions displayed in Fig. ^ 
embedded in a uniform perpendicular magnetic field B, 
which is a bipartite periodic structure with three sites 
per unit cell. As we shall see, the various characteristics 
of this system are similar to those discussed in Ref. || 
for two-dimensional tilings. Hereafter, we fix the total 
polarization to zero which is equivalent to consider two 
particles with opposite spin (f and J.). 

Let us consider the standard Hubbard hamiltonian : 



H = 
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where c\ a (resp. Cj iCT ) denotes the creation (resp. annihi- 
lation) operator of a fermion with spin a, = ct a t(J 
the density of spin a fermion on site i, and (...) stands 
for nearest neighbor pairs. Note that, since the particle 
considered here are fermions, the interaction term U is 
only efficient in the singlet sector where the orbital part 
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of the wave function is symmetric. When B = 0, the 
hopping term f y = 1 if i and j are nearest neighbors and 
otherwise. In the presence of a magnetic field O, tij 
is multiplied by a phase factor e llij involving the vector 
potential A : 

7 = / A.dl, (2) 

00 Ji 

where 0o — he/ e is the flux quantum. For convenience, 
we choose a gauge in which only one hopping term per 
unit cell is modified (see in Fig. [l]) . The whole spectrum 
only depends on the reduced flux / = 4>/4>q where <f> = 
Ba 2 /2 is the magnetic flux through an elementary square 
(a is the unit cell vector length). 




FIG. 1. Square chain under a uniform magnetic field. The 
magnetic flux per unit cell is denoted by <f) and 7 = 2n<f)/(j)o- 

Let us first analyze the one-particle problem. In this 
case, the translation invariance of the system along the 
chain direction allows one to straightforwardly compute 
the one-particle spectrum that consists of three bands : 

e a {k) = 2a y/l + cos(7/2) cos(fca) fce[0,27r/a], (3) 

where a = 0,±1 is the band index and 7 = 27r/. The 
weight of each band in the normalized density of states 
equals 1/3. The existence of a non dispersive band at 
e = is simply due to the bipartite character of the 
structure, its degeneracy being equal to the difference 
between the number of sites of each family. The most 
striking feature is that for / = 1/2 (half a flux quantum 
per unit square), the spectrum is made up of three non 
dispersive bands. As discussed in Ref. Q, this property 
leads to a complete lock-in of any wave packet spreading 
inside the so-called Aharonov-Bohm cages. One thus has 
a transition induced by the magnetic field. 

For {7 = 0, the two-particle spectrum is the addition 
of the one-particle spectra so that the eigenenergies are 
labelled by four quantum numbers : 

e a ^ ai (k h ki) = e QT (fc T ) +e ai (k l ), (4) 

where ct a = 0, ±1 (resp. k a ) is the band index (resp. the 
wave vector) of the spin a particle. Thus, for / = 1/2, 
the spectrum consists of five non dispersive bands corre- 
sponding to e = 0, ±2, ±4 and the space evolution of any 
two-particle wave function is confined in an Aharonov- 
Bohm cage that is merely the superposition of each one- 
particle cage. 



We now address the interacting case where U 0. The 
main question is whether or not the latter system remains 
frozen when the particles are interacting. In other words, 
can a (local) interaction term destroy the cages and au- 
torize any propagation ? In general, a two-particle prob- 
lem with on-site interaction in a D-dimcnsional structure 
can be viewed as a single particle one in a 2D-dimcnsional 
structure with a local potential in the hyperplane corre- 
sponding to a double occupancy of a site. Taking advan- 
tage of the translation invariance, this problem can then 
be mapped onto a (continuous) family of D-dimensional 
problems with a finite number of impurity sites which 
are often easier to solve O . The same approach could in 
principle be used here but given the very special nature of 
the non interacting system for f — 1/2, it is most easily 
carried out by using the minimally extended one-particle 
eigenstates displayed in Fig. ||. 
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FIG. 2. Eigenstates of the one-particle problem for 
/ = 1/2 (non normalized cage solutions). 

These states have non vanishing amplitude only on a 
finite number of sites and thus reflect the absence of prop- 
agation at / = 1/2. We denote them by \i,s a ), where 
i is the cell number centered on each 4-fold coordinated 
site and e a = 0, ±2. 

The two-particle state space is then given by the tensor 
product of the one-particle state : 

\i,a;i',a') = \i,e a )i ® \i',e a >)i, (5) 

for all i, i' , a, a' . It is worth stressing that since the \i, e a ) 
are confined eigenstates of the one-particle hamiltonian, 
most of these two-particle eigenstates are not affected by 
U. This local interaction term only acts on states for 
which the two particles have a non vanishing probability 
to be on the same site, i.e. such that \i — i'\ < 1. As 
a result, the number of states which are sensitive to U 
scales linearly with the total number of cells N whereas 
in the generic case (/ 7^ 1/2), it is proportional to A^ 2 . 
To proceed further, let us remark that the states sensitive 
to U are space-symmetric : 

\i, a; i', = — (\i, a; i' , a') + \i', a'; i, a}) , (6) 
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if i 7^ %' or a ^ a', and \i, a; i, a)s = \i, a; i, a)- Moreover, 
since the problem is invariant under a translation of the 
center of mass, it is convenient to build Bloch waves : 



JV-l 



\<p (a,a',K)) = -= > < 
v jv 



\n,a;n 7 a') s (7) 
<Pi(a,a/,K)) = 4= E e tKna \n,a;n + l,a') s , (8) 



n=0 
N-l 



n=0 

where K = 2tt j /Na(j = 0, TV— 1) is a wave vector lying in 
the first Brillouin zone. In the following, we shall always 
considerthe limit where iV tends to infinity. One then has 
to calculate the matrix element of the hamiltonian in each 
irreducible representation labelled by K, and restrict this 
analysis to the (15 x 15) subspace generated by \<po) and 
\<Pi). For any (a, a' , (3, (3' , K\, K m ) one has : 

{ip {a,a',K l )\H\<pi{P,P',K m ))=Q (9) 

{ipx{a t <J,Ki)\H \y x {(3,(3',K m )) = \ a , a ',p, > 6 t , m , (10) 

where \a,a',p,0' is a _fT; !m -independent scalar and where 
3i,m is the usual Kronecker symbol. This implies that the 
eigenvalues of the (9 x 9) subspace generated by the \ipi) 
are non dispersive. These eigenvalues are given by the 
five [/-dependent roots of the characteristic polynomial : 

P(e, U) = e 5 - Ue 4 - 20 e 3 + !6Us 2 + 64e - 24Z7, (11) 

and four [/-independent values ±2, (two-fold degene- 
rated) resulting of additional symmetry between a and 
a' . The non dispersive part of the spectrum 113 is shown 
in Fig. I for U > 0. 
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FIG. 3. Non dispersive part of the two-particle spectrum 
at / = 1/2 versus the interaction. 

We emphasize that since the spectrum of the Hubbard 
hamiltonian (in a bipartite, structure) is odd under the 
transformation {/—>—{/ M, we can restrict our analysis 



to the repulsive case. In the large U limit, the five roots of 
P tends toward ± v8± 2\/l0, U ; this latter value simply 
corresponding to a situation where the two particles are 
localized on the same site (anti-bonding state). 

A much more interesting component of the spectrum 
is provided in the (6 x 6) subspace «LB generated by the 
Ivo)- I n ^ n is subspace, the eigenvalues are given by the 
roots of the following characteristic polynomial : 

Q(e, U, K) = e 6 - 2Ue 5 + (U 2 - 20) e 4 + 28Ue 3 + 

8 (8 - U 2 ) e 2 - 4U(U - 3 cos(Ka))e + 
4U 2 (2 + cos(Ka)). (12) 

Contrary to the previous case, these eigenvalues obvi- 
ously depend on K, and the associated non degenerated 
eigenstates are extended (Bloch-like). This dispersive 
component of the spectrum is displayed in Fig. 
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FIG. 4. Dispersive part of the two-particle spectrum at 
/ = 1/2 versus the interaction. 

In the large U limit, the asymptotic eigen- 
values are given by U (twice degenerated) and 



± y 4 ± 2 -y/2 — cos(Ka) , and there still remains a in- 
dependent component in the spectrum. 

The physical consequences of this dispersion are im- 
portant. Indeed, let us consider a generic two-particle 
initiaLstate having a non zero overlap with one of the 
\ipo) t3. The emergence of dispersive states for U 7^ 
indicates that it is now possible for this wave packet to 
spread over the whole system whereas it was completely 
trapped inside the Aharonov-Bohm cage in the non in- 
teracting case. Moreover, since the dispersive eigenstates 
are extended, the propagation is ballistic. 

For other values of the reduced flux, the full solution of 
the two-particle problem cannot be cast in a simple ana- 
lytic form, although it is possible to reduce it into a scat- 
tering problem for one particle moving on a chain in the 
presence of three static impurities. Nevertheless, let us 
analyze the neighborhood of the half-flux parametrized 
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by Sf = |/ - 1/2|. For / ^ 1/2, all the single-particle 
eigenstates become extended, except those corresponding 
to the flat band at e = 0. The non dispersive two-particle 
states which are insensitive to U at / = 1/2 evolve into 
several two-particle continua, whose band width scales 
as Sf for small Sf. The non dispersive states which are 
sensitive to U become dispersive with a band width scal- 
ing as (Sf) 2 /U (see inset in Fig. ||). The corresponding 
wave functions still exhibit a binding of the two parti- 
cles. Finally, as displayed in Fig. |j| the dispersive states 
at / = 1/2 remain dispersive. All these dispersive bound 
states evolve smoothly as a function of Sf until their en- 
ergies merge in the two-particle continuum, which occurs 
for Sf ~ U at small U. As Sf further increases, the total 
number of bound states gradually decreases. 
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FIG. 5. Low-energy spectrum for U = 0.01 as a function 
of the reduced flux. 

Since, in this study, the effect of the interaction term 
between the two particles is clearly to induce a delocaliza- 
tion process, one can expect that for finite density of par- 
ticles, a subtle correlated conducting state will emerge. 
Note however that in this system, the energy of two parti- 
cles is minimal when they do not interact, either because 
they are far apart or because their orbital wave function 
is antisymmetric. We therefore conjecture that the many 
body ground state will remain localized up to a filling fac- 
tor equal to 1/3. At this point, the lowest flat band is 
completely filled with a fully polarized electron sea. The 
next electron will have an opposite spin and is thus likely 
to delocalize along the chain. This problem obviously 
deserves futher investigations. Let us also remark that, 
here, we consider a quasi-one-dimensional model to sim- 
plify the calculations, but a similar physics is expected for 
other tight-binding models which exhibit single-particle 
confinement inside Aharonov-Bohm cages. Nevertheless, 
in more general systems, we do not expect to see non 
dispersive states sensitive to U. In this context, an inter- 
esting open question is whether it is possible to find situa- 



tions where the ground state of the two-particle spectrum 
is dispersive, even for repulsive interaction. 

How could this interaction induced delocalization be 
observed experimentally? The two-particle system might 
be accesible in a Josephson junction or a quantum dot 
array in the Coulomb blockade regime. Such experimen- 
tal situations shall be described by a model with on- 
site disorder, but, very likely, some two-particle states 
would still exhibit a much larger localization length than 
single particle ones. Ongoing experiments on ballistic 
semi-conducting networks with special two-dimensional 
geometry fe3 may also manifest some interaction effects 
on this Aharonov-Bohm localization. In these structures, 
the interaction strength can be varied by changing the 
electronic density or by polarizing the system with a 
tilted magnetic field. Clearly, a better understanding of 
the many electron case needs to be achieved. 

We would like to thank CI. Aslangul, B. Delamotte, 
D. Mouhanna and M. Tissier for fruitful discussions. 
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The dimension of this subspace is different from the one 
generated by the since, for any (a,a',K) one has : 
\ip {a,a',K)) = \ip (a ,a,K)). 

Such a situation occurs as soon as the two one-particle 
wave functions overlap each other. For example, one can 
considers a initial configuration where the two particles 
are localized on the same site. 
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